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In polyatomic molecules with II electronic ground state the ro-vibrational spectrum can be 
strongly modified by the Renner- Teller effect. The linear form of C3H molecule has particularly 
strong Renner- Teller interaction and a very low lying vibronic level, which corresponds to the 
excited bending vibrational mode. This leads to the increased sensitivities of the microwave and 
submillimeter transition frequencies to the possible variation of the fine structure constant a and 
electron to proton mass ratio [i. 



PACS numbers: 06.20.Jr, 06. 30. Ft, 33.20.Bx 
I. INTRODUCTION 

At present it is generally recognized that microwave 
and submillimeter molecular spectra from the inter- 
stellar medium provide us with a very sensitive tool 
to study possible variation of the fundamental con- 
stants a — e 2 /Tic and /i = m e /m p on a large space- 
time scale. It was shown that certain types of transi- 
tions are particularly sensitive to such variations. The 
mixed tunneling-rotational transitions in such molecules 
as H 3 0+, H 2 2 , CH3OH, and CH 3 NH 2 can be very sen- 
sitive to /i- variation [THE]. Recently the submillimeter 
spectra of methanol have been used to place very strin- 
gent limits on /i-variation on the cosmological timescale 
[5J[7]- On the other hand, the A-doublet transitions in 
such diatomic radicals, as OH and CH are very sensitive 
to variation of both constants [BTUOj. The 18 cm OH line 
was observed at high redshifts, which allowed to constrain 
time variation of both constants [11 . In that work the 
21 cm hyperfine hydrogen line was used as a reference. 
This constrain can be further improved if more than one 
A-doublet transition in OH, or CH is detected. 

A rather general way to look for the enhanced sensitiv- 
ity to variation of the fundamental constants is to search 
for the accidental degeneracy of the levels of different 
nature. This approach works for very different systems 
from nuclei, to atoms, and molecules (see, for example, 
the review [T2] and [IB"]). 

In this paper we want to draw attention to the mi- 
crowave and submillimeter spectra of the linear poly- 
atomic radicals with nonzero electronic angular momen- 
tum. First, these molecules have iv-doublets, which are 
analogous to the A-doublets in diatomics. Second, the 
Renner- Teller interaction can lead to the anomalously 
low lying vibronic levels and cause enhanced sensitivities 
of the mixed ro-vibronic transitions. Finally, there are 
many linear polyatomic molecules, which are detected in 
the interstellar medium. In this context one of the most 
interesting species is the linear C3H molecule, where the 
excited vibronic level lies less than 30 cm" 1 above 
the ground level TI1/2 and where several mixed transi- 



tions were recently measured in the molecular beam ex- 
periment [H] . Because of that we focus on this molecule 
here leaving other similar molecules for a separate dis- 
cussion. 

The interstellar carbon-chain radicals of hydrocarbon 
series C„H (n — 2 — 6) exist in linear and cyclic isomeric 
forms. Both forms are observed in the millimeter-wave 
range toward dark and translucent molecular clouds and 
circumstellar envelopes of carbon-rich stars [T5H2T] . A 
typical abundance of the linear radical I-C3H, — the 
simplest odd carbon chain radical under consideration 
in the present study, — is ~ 10 -9 relative to hydro- 
gen. The cyclic-to-linear abundance ratio [c— C 3 H]/[/- 
C3H] ~ 5—10 in dark clouds [23l[25], but decreases to ~ 3 
in translucent clouds |23j . and becomes less than unity 
around carbon stars |27] . The cyclic and linear isomers of 
C3H have also been detected in a star-forming region [23] 
and in two extragalactic sources: the star-burst galaxy 
NGC 253 [25] and the spiral galaxy located in front of the 
quasar PKS 1830-211 at the redshift z = 0.89 [3D]. Thus, 
I-C3H lines have been detected in regions with kinetic 
temperature ranging from T^m ~ 10 K (dark clouds) 
to several hundred Kelvin (circumstellar envelopes, star- 
forming regions). The observed line intensities are less 
or about 0.1 K. 

The lines observed from the interstellar medium (ISM) 
are Doppler broadened, so the linewidth r « Trj = u> , 
where AV is the velocity distribution for the ISM and c 
is speed of light. This means that - « characterizes 
ISM and to a first approximation is independent on the 
frequency of the transition u>. Because of that for the 
astrophysical search of the possible variation of the fun- 
damental constants it is crucial to find transitions with 
high dimensionless sensitivity coefficients defined as: 

Suj „ See 811 

— = Q a — + Qn—- (l) 
uj a 11 

In the optical waveband these sensitivity coefficients are 
typically of the order of 10 -2 . In the microwave wave- 
band they are typically of the order of unity, but can 
be much bigger. Below we will calculate Q factors for 
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some microwave and submillimeter transitions of the l- 
C3H molecule and show that they can reach the order 
10 3 . As usual such enhanced sensitivities take place for 
the low frequency transitions between quasi degenerate 
levels of different nature. 



H so = A cS AE , A cS = AA eS /A . 



(3) 



Let us define the model more accurately. Following 
1341 we write Hamiltonian as: 



II. RENNER-TELLER EFFECT 

In this section we briefly recall the theory of the 
Renner- Teller effect in polyatomic linear molecules [3IJ 
[32"] . The total molecular angular momentum J consists 
of several contributions. We have the overall rotation 
of the nuclei R. Then we have the vibrational angu- 
lar momentum G associated with the twofold degener- 
ate bending vibration mode(s) and the electronic angu- 
lar momentum L. Momentum R is perpendicular to the 
molecular axis £; two others have £ projections I and A. 
We define momentum TV = R + G + L and its projec- 
tion (Nq) = K = I + A. Finally, we add electronic spin: 

j = at + s, ( j c > = n. 

Suppose we have II electronic state |A = ±1) and v = 
1 vibrational state of a bending mode \l = ±1). All 
together there are 4 states |A = ±1)|^ = ±1). We can 
rewrite them as one doublet A state \K = ±2) and states 
S + and In adiabatic approximation all four states 
are degenerate. Renner [3T] showed that the states with 
the same quantum number K = I + A strongly interact, 
so S + and £~ states repel from each other, while A 
doublet in the first approximation remains unperturbed. 
We are particularly interested in the case when one of 
the £ levels is pushed close to the ground state v = 0. 
This is what takes place in Z-C 3 H molecule [T4"HT9"1 133] . 

Consider linear polyatomic molecule with unpaired 
electron in the state in the molecular frame £, 77, £. 
Obviously, the bending energy is different for bendings in 
and in 7?£ planes: V± — \k±x 2 (here \ is the supple- 
ment to the bond angle) . That means that the electronic 
energy depends on the angle 4> between the electron and 
nuclear planes: 



H — H„ + T„ 



AL C S, 



H' = V' cos 2(j> , 



(2) 



where 2V' = V+ — V- = k'x 2 ■ There is no reason for V' to 
be small, so k' ~ k± ~ 1 a.u. and to a first approximation 
k' does not depend on a and /1. 

As long as interaction ^ depends on the relative an- 
gle between electron and vibrational rotation it changes 
angular quantum numbers as follows: AA = —Al = ±2 
and A AT = 0. This is exactly what is necessary to pro- 
duce splitting between S + and £~ states with v = 1 as 
discussed above. 

Interaction |2]) also mixes different vibrational levels 
with Aw = ±2, ±4, .... Thus, we have, for example, 
the nonzero matrix element (ME) (0, 0, 1, 1\H'\2, 2, -I , I) 
between states \v, I, A, K). Such mixings reduce effective 
value of the quantum number A and, therefore, reduce 
the spin-orbital splitting between II4/2 and TI3/2 states 



(4) 



Here "electronic" part H e includes all degrees of freedom 
except for the bending vibrational mode and spin. For 
I-C3H there are two bending modes, but for simplicity 
we include second bending mode in H e too. Electronic 
MEs in the |A) basis have the form: 



(±1 \HJ ± 1) = 



(±MHe\Tl) 



k' 



V- _ k 2 
~2 ~ 2 X 

X 2 exp (=p2z0) 



(5a) 
(5b) 



Here x an d <j> are vibrational coordinates for the bending 
mode. Kinetic energy in these coordinates has the form: 



T 



f 



if 



1 d 



1 d 2 



2MR 2 \d X 2 xdx X 2 d<j> 2 



(6) 



We can use the basis set of 2D harmonic functions in 
polar coordinates p = X-R and <j> for the mass M and 
force constant k: 



ipv,l (A 0) = Rv,i (p) —7= exp (il(j>) 



2ir 



(7) 



It is important that radial functions are orthogonal only 
for the same I: 



(Rv',i\Rv,i) = S V ' tV . 



(8) 



This allows for the nonzero MEs between states with dif- 
ferent quantum number I. Averaging operator Q over 
vibrational functions we get: 

(v',l'\H e + T v \vJ) = [u v (v + l)+AAS c ]6 v r tV 8 V j 

+ \{Rv>v\k' x 2 \R v i) exp {T2i<j>)8 v> i ± 2 . (9) 

The exponent here ensures the selection rule A' = A =p 2 
for the quantum number A when we calculate MEs for 
the rotating molecule. 

We solved eigenvalue problem for Hamiltonian Q us- 
ing the basis set of the 2D-harmonic oscillator. Matrix 
elements were formed according to Eq. (|9| . As discussed 
above we neglected one of the bending modes leaving only 
the one that produces K — level close to the ground 
state doublet K = 1, fi = 1/2,3/2. Our model Hamil- 
tonian has only 3 parameters, namely u v , A, and the 
dimensionless Renner- Teller parameter £: k' = £k. In 
Ref. [31] the following values were obtained: 



u> v 



589 cm" 1 , A = 29 cm" 1 , £ = 0.883 . (10) 
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We fixed the values for u„ and A and varied the Renner- 
Teller parameter £ to fit five lowest levels for the given 
bending mode: n 1/2 , n 3/2 , S 1/2 , A 3/2 , and A 5/2 . The 
optimal value appeared to be £ = 0.788. The difference 
with ( 10 ) is probably due to the neglect of the anhar- 



monic corrections and second bending mode. 

Our results are presented in Table |T] The first two 
columns give nominal vibrational quantum number v and 
its actual average value. We see that Renner-Teller term 
in Q strongly mixes vibrational states. This mixing 
also affects (A) and decreases spin-orbital splittings as 
explained by Eq. 



TABLE I: Low lying energy levels for the bending mode w v = 
589 cm -1 and their sensitivities q a and q^ to the variation of 
a and fj, respectively. All values are in cm -1 . 



V nom 


(v) 


K 


Q, 


(A) 


E A 


34 


9m 


q a 





1.22 


1 


0.5 


0.50 


367.9 0.0 


0.0 


187.8 


-14.6 





1.35 


1 


1.5 


0.46 


381.9 13.9 


14.0 


187.8 


13.3 


1 


2.32 





0.5 


-0.01 


394.2 26.3 


27.0 


197.3 


-0.4 


1 


3.57 


2 


1.5 


0.21 


597.7 229.7 


226.0 


300.3 


-6.1 


1 


3.65 


2 


2.5 


0.19 


603.5 235.5 


232.0 


300.3 


5.5 



The last two columns in Table |T] give sensitivity coeffi- 
cients q a and g M in cm' 1 : 



SE = q a - 



1n- 



To get them we assumed that parameters ( 10 1 scale in 
a following way: ui v ~ y 1 ^ 2 , A ~ a 2 , and £ does not 
depend on a and y,. The dimensionless sensitivity coeffi- 
cients for the transitions cj^k = E^—Ei can be found 
as: 

Qi,k = (ik - ?t)M,* • 

In Table [TT] these coefficients are calculated for the same 
set of parameters as in Table [IJ and for the slightly differ- 
ent parameters which better fit experimental frequencies 
from Ref. [13] • We see that Q- factors are practically the 
same for both sets. 



TABLE II: Q-factors for the transitions between states from 
Table [IJ and for parameters A e ff and AEy.ii defined by (|3| and 
|TTJ respectively. Frequencies are in cm -1 . 



K Q 



0.5 1 
1.5 
0.5 2 
1.5 2 

A C ff 
A£ En 



k' ty 

1.5 
0.5 
1.5 
2.5 



Fit to Ref. |34j 



Fit to Ref. 

W Qa 
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13.9 
12.4 
203.5 
5.8 
13.9 
19.4 



0.00 
0.78 
0.51 
0.00 
0.00 
0.50 



2.00 14.4 

-1.11 13.3 

-0.03 204.4 

2.00 6.0 

2.00 14.4 

0.00 20.5 



0.00 
0.77 
0.51 
0.00 
0.00 
0.50 



2.00 
-1.07 
-0.03 
2.00 
2.00 
0.00 



is significantly reduced by the Renner-Teller mixing: the 
fine-structure parameter is 29 cm -1 and the splitting be- 
tween IIj/2 and n 3 / 2 is only 13.9 cm" 1 . According to 
% the mixing reduces the splitting. However, this ef- 
fect depends on the dimensionless Renner-Teller param- 
eter £ and does not depend on /i and a. Consequently, 
the effective parameter A e s depends on the fundamental 
constants in the same way as initial parameter A. 

For the high frequency transition £1/2 — ► A3/2, where 
spin-orbital energy can be neglected, we get Q M = 0.5 and 
Q a = 0. These results are expected, because our model 
has only two dimensional parameters: vibrational fre- 
quency, which is proportional to /1 1 / 2 and the fine struc- 
ture parameter A, which scales as a 2 . Even though our 
vibrational spectrum is far from that of a simple har- 
monic oscillator, the non-diagonal MEs ^ of the Hamil- 
tonian Q still scale as y}l 2 . Therefore, if we neglect 
spin-orbital splittings, we get = 1/2 for all transi- 
tions. The only transition in Table [TT] where spin-orbital 
energy and vibrational energy are close to each other is 
the II3/2 — > £1/2 transition. Resultant frequency is 
roughly half of the vibrational energy difference between 
II and £ states. This leads to s» 1 and Q a ps — 1. 



The following analysis in Sec. Ill will be based on the 
effective Hamiltonian for the rotating molecule. The lat- 
ter includes only two parameters from this section: the 
effective fine-structure parameter A e g and the energy dif- 
ference between E and II states, 



AE- 



sn 



E(Tl 1/2 ) + E(Tl 3/2 ) 



(11) 



Numerical values for these parameter will be obtained 
from the fit to experimental transition frequencies. Here 
we only need to determine the dependence of these pa- 
rameters on fundamental constants. Table iHl shows that 
A e s ~ a 2 and A£^n ~ /J 1 / 2 . Once again, this is because 
the Renner-Teller mixing depends on the dimensionless 
parameter £ and does not depend on a and y. 



III. EFFECTIVE HAMILTONIAN FOR 
ROTATING MOLECULE 

In this section we mostly follow Ref. [T3] . However, we 
prefer to use the basis set for the Hund's case "a" as we 
did before in [10, 35 . We define the effective Hamiltonian 
for the subspace of the three lowest vibronic states IIi/ 2 , 
II3/2 j and SJ"/ 2 . We neglect hyperfine interaction and 
some minor centrifugal corrections included in [14] . 

The basis ro-vibronic states for the Hund's case "a" 
have the form: 

\v, I, A, (K), S, S, J, n, M) = \v, I) |A) \S, E) | J, n, M) . 



For the two fine structure transitions, 
and Act/9 — ► A 5 / 2 , we get sensitivities Q M 



-> n 3/2 

and 



3/2 " 

2. This may seem strange as the fine structure 



Here the quantum number K does not appear explicitly, 
being defined as K = I + A. From these basic states we 
form parity states as described in [36] : 
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|n va ) = 



n 3/2 > = 



|0,0,1,(1), i,-±,J,l,M,p) = i|0,0}(|l}|i,-i)|J ! i,M)+Xp|-l}||,|)|J,-|,M}) , 



10,0,1,(1), §, |, J, §,M,p) = i|0,0)(|l>|i |>|J, |,M> + X p| - l)|i,-i>|J,-|,M)) , 



(12) 
(13) 



]+ /2 ) = |1, 1, 1, (0), 1, 1 J, l,M,p) = i(|l, 1)| - 1) + |1, -1)|1)) (||, |)| J, |,M) + x p |A,-i)|J,-§,M}) , (14) 



j-s 



where the parity dependent phase is \p — ( — 1) 

We can write rotational energy by adding vibrational 
angular momentum G to the usual expression: 

H rot = B(J -G-L-S) 2 = B[J(J + 1) - n 2 } (15a) 

— 2B [JqGq + JqL q + J q S q 

q=±l 

+ G q L_ q + G q S- q + L q S- q ] (15b) 
[GqG-q + L q L_ q + S q S- q ] . (15c) 

q=±l 



Here we use the recipe from that in the molecular 
frame all scalar products involving total angular momen- 
tum J are written as J q X q rather than (— l) q J q X- q . The 



last line (15c) can be skipped as it gives a constant inde- 



pendent of J, il, and p. The terms in (15b) linear in L q 
turn to zero in the subspace A = ±1. We are left with 
the following operator for the rotational energy: 

H Tot = b[j(j + 1) - n 2 } - d[j(j + 1) - n 2 } 2 

- 2B [JqGq + JqS q + G q S- q ] , (16) 

q=±l 

where we added standard centrifugal correction to the 
main diagonal term. 

It is straightforward to calculate MEs of this operator 



on the states ([12]) - ( 14 ). The term J q S q does not change 
quantum number I and can not mix E and II states. The 
nonzero matrix elements are: 



3/2 1 



- 2BJ a S, 



|ni /2 ) = -b^(j-±)(j + 1), (17) 



E+ \-2BJ g S q \X+ ) = -B Xp {J+\ 



(18) 



The operator J q G q changes quantum number / by one 
and mixes E and IT states: 



J l/2l 



-2BJ q G q \Tl l/2 ) =-(i Xp {J- 



(E+ a |-2BJ,G,|n 3/2 ) = /y (J - §)(J + f) ■ 
where /3 is defined as 

p = B(l = l\Gx\l = 0) . 



(19) 
(20) 

(21) 



This ME can not be calculated within this formalism and 
is included as an independent parameter of the effective 



Hamiltonian (see also Sec. IV). Finally, the term G q S 



mixes E and II states, but can not change the quantum 
number f2: 



(T,^, 2 \ — 2BG q S- 



|n 1/2 ) = -/?. 



(22) 



In addition to the rotational energy the effective 
Hamiltonian must include spin-orbit interaction ([3]), the 
energy splitting between E and n states Ai^sn and spin- 
rotation interaction. Following [35] we write the latter as: 

1 (NS)= 1 (J-S)S 

= 7 |f!E +J2 JqS q - S{S + 1)J . (23) 

The nontrivial part of this interaction is now reduced to 



the MEs (17) and (18) 



Equations (19) - (22) show that the Coriolis terms 



involving vibrational angular momentum G lead to the 
if-doubling via interaction between n and E states. In 
contrast to the terms involving electronic angular mo- 
mentum L here we do not need mixing with excited elec- 
tronic states. Still, because of the relative smallness of 



the parameter f3 in (21 ), these latter terms can not be ne- 



glected. They have exactly the same form as for diatomic 
molecules and are defined in Ref. |36j . 

Transition amplitudes between spin-rotational states 
of the l-G^H molecule are expressed through MEs of the 
dipole moment operator D on the basic states ( 12 ) - 
(14). Generally speaking there are both diagonal and 



nondiagonal MEs in vibrational quantum numbers v, I. 
Let us estimate them using atomic units (h = m e = 
|e|=l). 

In the molecular frame the diagonal ME is reduced to 
the dipole moment of the molecule along the molecular 
axis (v,l\D^\v,l) « D. If we assume that the charge 
of the hydrogen atom in the molecule is q, then D ~ 
2qRo ~ 4q, where Rq is the bond length. Comparing 
this estimate with calculated value D = 1.40 [37] we get 
q = 0.35. Now we can estimate the nondiagonal ME: 
(0,0|Di|l,-l) ~ q£ ~ q/y/Mu^ ~ gM" 1 / 4 ~ Q.lq ~ 
0.04, where £ is the amplitude of the vibration and M ~ 
10 4 is the reduced mass for this vibration mode. We 
conclude that nondiagonal MEs are much smaller than 
diagonal, so we will neglect them. 

In this approximation we get the following expressions 
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for the reduced MEs on the basis states (12) - (14) 



(Xn,J',p'\\D\\Y n ,J,p)=6 x , Y (-l) J '- il 

J' 1 — p'p 

-non 



xV(2J' + l)(2J+l)( J fi I a) 



D, (24) 



where X and Y denote either II, or £ state. Below we use 
these expressions and theoretical value D = 1.40 a.u. [37] 
to estimate reduced MEs for the microwave transitions 
in Z-C3H. The Einstein coefficients A for these transitions 
can be found as 1381: 



Ai-yj — 



Hm\D\\j)\ 2 ^ 



(25) 



3kc 2 Ji + 1 
where reduced ME is in a.u. and do is the Bohr radius. 



IV. SCALING OF THE PARAMETERS OF THE 
EFFECTIVE HAMILTONIAN WITH a AND fx 

Effective Hamiltonian described in Sec. IIIII is essen- 
tially equivalent to the one used in Ref. [TJ]. We in- 
cluded centrifugal corrections to most of the terms using 
the same definitions as in [14]. For the hyperffne struc- 
ture we used usual parameters a, 6, c, and d. Note that 
in [T3] the constant bp = b + c/3 was used instead of b. 

In this section we discuss how the parameters of the 
effective Hamiltonian depend on the constants a and fi 
(see Table III). The scaling of the two largest parameters, 



A-Esn ~ fi 1/z and A e g ~ a , has been already discussed 
in Sec. [n] The rotational constants By and Bn linearly 
depend on /1. The spin-rotational interaction (23) ap- 
pears from the second order cross term in Coriolis and 
spin-orbit interactions, therefore 7 ~ a 2 /1. For n states 
there are two additional terms of the spin-rotational in- 
teraction with parameters p and q. The first of them has 
the same scaling, as 7, i.e. p ~ a 2 /i. The second term is 
quadratic in Coriolis interaction, so g ~ /1 2 . These scal- 
ings are obvious from the expressions on pp. 362 and 531 

of [36]. 

Let us now discuss the parameter (3 defined by (21 ). It 



is proportional to the nondiagonal ME (l|Gi|0). Accord- 
ing to Eq. (13) in [32], the perpendicular component G\ 
simultaneously depends on the vibrational coordinates of 
the bending (vb) and stretching (v s ) modes. In the har- 
monic approximation it has nonzero MEs only between 
different stretching vibrational states, i.e.: 

(v b = l,v a = l\Gx\v b = Q,v a = 0) ^0. 



The ME in (21) is diagonal in stretching quantum num- 
ber v s . It is nonzero due to the anharmonic corrections 
which mix vibrational modes. Such corrections appear in 
the first order of the adiabatic perturbation theory and 
are proportional to the adiabatic expansion parameter 
/1 1 / 4 . Thus we can expect that /3 ~ OAB. This estimate 
agrees well with the numerical value obtained in Sec. [V] 
We conclude that p ~ /// 4 B - /1 5 / 4 . 



TABLE III: Parameters of the effective rotational Hamilto- 
nian and their scaling with a and p. 



Par am. 


This work 


Ref. [Til 


Units Scaling 




609.9811 


609.9742 


GHz 


aV /2 


Be 


11.2124327 


11.2126703 


GHz 


aV 


D s 


4.548 


4.867 


kHz 


«V 


7s 


-35.800 


-35.525 


MHz 


aV 


7s, D 


18.04 


0.549 


kHz 


o-V 


&£ 


-6.3 


-6.29 


MHz 


aV 


CE 


31.8 


27.17 


MHz 


aV 


A cB 


432.7762 


432.7898 


GHz 


aV 


B n 


11.1892055 


11.1891033 


GHz 


aV 


D n 


5.356 


5.2340 


kHz 


aV 


7n 


—48.652 


—48.075 


MHz 


aV 


7n,D 


21.670 


0.000 


kHz 


aV 


P 


—6.9021 


—7.0681 


MHz 


«V 


1>D 


—1.595 


0.504 


kHz 


aV 


q 


-12.8556 


-12.9922 


MHz 


aV 




-0.443 


-0.1432 


kHz 


aV 


p 


1.2586 


1.2342 


GHz 


aV 5/4 




-28.3 


-19.2 


kHz 


aV 9/4 


an 


12.43 


12.32 


MHz 


aV 


6n 


-22.57 


-23.04 


MHz 


aV 


cn 


27.56 


28.07 


MHz 


aV 


dn 


16.21 


16.26 


MHz 


a^ 1 



Our effective Hamiltonian includes centrifugal correc- 
tions {D, 7/), /3d, etc.) to the most important terms. 
We assume that such corrections have the same a depen- 
dence as the respective main term and an extra power in 
their /i dependence. The magnetic hyperhne constants 
scale as the product of the nuclear and electronic mag- 
netic moments, i.e. as a 2 ^,. 

All scalings discussed above are approximate. There 
are relativistic corrections to all parameters, which mod- 
ify their a-dependence. These corrections are of the or- 
der of (aZ) 2 ~ 0.2%. The fj, dependence of parameters is 
changed by non-adiabatic corrections. To illustrate this 
point let us consider the rotational constants B. To a 
first approximation the small difference between By, and 
-Bn can be related to the vibrational corrections to the 
adiabatic value of the rotational constant Bq. 

We can use the data from Table Q] and Table IIIII to es- 
timate vibrational correction to the rotational constant: 



B v =B - a(v + 1) , 



B t 



Bn — By, 

vy - vn 



-22.8 MHz. 



(26) 

11137.1MHz, (27) 
(28) 



If we assume that Bq scales as \x and a scales as /z 3//2 [S] , 
we get following scalings of the rotational constants from 
Table Ell 



By- i-i 



1.010 



Bn ~ A« 



1.007 



(29) 
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Note that we neglected other vibrational degrees of free- 
dom, so actual corrections can be somewhat bigger. We 
conclude that we know the scalings of the main param- 
eters from Table |III| roughly to a percent accuracy. Fur- 
ther improvement of this accuracy requires extensive ah 
initio calculations. 



V. NUMERICAL RESULTS FOR ROTATING 
MOLECULE 



1400 - 




'%/2 



FIG. 1: Spin-rotational levels of the three lowest vibronic 
states of the molecule Z-C3H. if-doubling is shown schemati- 
cally. The levels are labeled by the quantum numbers J and 
p. The hyperfine structure for the two lowest if-doublets is 
shown on the inset. Due to a strong Renner- Teller effect the 
component 2 E + of the excited bending state 1/4 (CCH bend- 
ing) is shifted towards lower energies, ~ 29 cm -1 above the 
zero-level of the ground state 2 iI 1/ / 2 - 



TABLE IV: Frequencies (MHz), Q-factors and reduced MEs 
(a.u.) of some rotational transitions for Ifi/2, If3/2, and 
states. 





J-> J 4 


- f transitions for IL/2 state 




J p 


Jp 




Q ex 






ll^ll 2 


5 + 
2 


1 

2 


32627.84 


0.06(0) 





97(1) 


1.86 


3 
2 


1 + 


32662.10 


0.06(0) 





97(1) 


1.86 


5 
2 


1 + 


54405.75 


0.06(0) 





97(1) 


3.33 


1 + 


3 
2 


54428.34 


0.06(0) 





97(1) 


3.33 


1 + 


5 
2 


76199.10 1 " 


0.06(0) 





97(1) 


4.72 


7 

2 


2 ^ 


76204.62+ 


0.06(0) 





97(1) 


4.72 


35 
2 


33 , 

2 


383435.12* 


0.02(0) 





99(1) 


18.51 


35 , 

2 "r 


33 
2 


383942.45* 


0.02(0) 





99(1) 


18.72 


47 
2 


45, 

2 ~ 


516312.64* 


0.01(0) 





99(1) 


22.39 


49 
2 


47, 

2 ~ 


539280.62* 


0.01(0) 





99(1) 


23.43 




J —> J + 1 transitions for II3/2 state 




J'P' 


Jp 




Qa 






ll-DII 2 


5 
2 


2 ^ 


57437.17 


-0.05(1) 


1 


03(1) 


2.22 


1 + 


3 
2 


57453.77 


-0.06(0) 


1 


03(1) 


2.22 


f + 


5 
2 


80388.41 


-0.05(0) 


1 


03(1) 


3.93 


7 
2 


2 T 


80421.07 


-0.05(0) 


1 


03(1) 


3.93 


47 
2 


45 , 
2 ^ 


532658.83* 


-0.06(0) 





98(1) 


17.54 


49 
2 


47, 

2 ~ 


556392.87 


-0.01(0) 


1 


00(1) 


23.16 


49 , 
2 ^ 


47 
2 


552385.85* 


-0.12(0) 





94(1) 


13.61 


51 
2 


49, 
2 ^ 


599557.31 


0.09(1) 


1 


00(1) 


14.29 


51 , 

2 ~ 


49 
2 


578834.56 


-0.01(0) 


1 


00(1) 


23.77 


53 
2 


— + 

2 ~ 


601263.86 


-0.01(0) 


1 


00(1) 


24.38 


TV — > N + 1 transitions for E+ /2 


state 




N' J'p' 


N Jp 




Qot 






ll^ll 2 


1|- 


0| + 


22468.22 


0.00(0) 


1 


00(1) 


0.93 


1|- 


0i + 


22420.66 


0.00(0) 


1 


00(1) 


1.87 


2| + 


1|- 


44888.43 


0.00(0) 


1 


00(1) 


1.87 


2| + 


1 | - 


44857.31 


0.00(0) 


1 


00(1) 


3.36 


23 f - 


22 f + 


517808.87* 


0.05(0) 


1 


02(1) 


26.02 


24 f + 


23 f - 


538178.20 


0.00(0) 


1 


00(1) 


32.57 


24 f + 


23 f - 


542975.34 


0.11(0) 


1 


07(2) 


19.58 


25 f - 


24 f + 


560575.62 


0.00(0) 


1 


00(1) 


33.93 


25 f - 


24 f + 


540680.86 


-0.11(1) 


1 


01(2) 


18.11 


26 f + 


25f- 


586696.44 


-0.07(0) 





94(1) 


23.27 



Transitions detected at the redshift z = 0.89 in Ref. [30] . 
* Transitions observed in Ref. [T¥|. 



Our effective Hamiltonian has 22 parameters listed in 
Table |III| including 6 parameters for the hyperfine struc- 
ture. The 16 non-hyperfine parameters were fitted us- 
ing simplex method to the 44 experimentally observed 
transitions from Ref. [14] and to 12 experimental fre- 
quencies listed in the NIST database [33]. We added 
8 theoretically predicted transitions for lower rotational 
quantum numbers from the same database to be sure we 
adequately reproduce this part of the spectrum. 

In our fit the rms deviation for 64 fitted transition is 
0.23 MHz with maximum deviation 0.52 MHz. This ac- 
curacy is lower than typical accuracy of the similar fits 
in the literature, but is absolutely sufficient for our pur- 
poses. Our main goal here is to calculate sensitivity co- 
efficients for different transitions to the variation of the 



fundamental constants. Though the sufficiently complex 
effective Hamiltonians allow for very accurate predictions 
of the transition frequencies, the accuracy they can pro- 
vide for the sensitivity coefficients is limited by the un- 
certainty in the dependence of the used parameters on 
the fundamental constants (see Sec. IV and [S]). 

To fit the hyperfine structure parameters we used 30 
lines from the Ref. [19] and 12 i^-doublet transitions 
from [ID] . The hyperfine structure is mostly too small to 
change the values of the sensitivity coefficients. This in 
not true only for several -ftT-doublet transitions with fre- 
quencies < 100 MHz, comparable to the hyperfine split- 
tings. We used the scalings from Table |III| to calculate 
the shifts of the spin-rotational levels due to the change 
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of the constants a and [i by ±0.1%. After that we found 
dimensionless sensitivities Q a and for the transitions 
described by the effective Hamiltonian. 



TABLE V: Frequencies (MHz), Q-factors and reduced MEs 
(a.u.) for A'-doubling transitions in II1/2 state. 



JF'p',Fp 


u 


Qa. 




WW 




0- 


52.37 


0.66(2) 


1-7(2) 


0.333 


|0+, 


1- 


39.12 


0.20(2) 


1.9(2) 


0.333 




1- 


34.93 


-0.02(2) 


2.0(2) 


0.667 


|1- 


1+ 


85.55 


0.65(2) 


1.7(1) 


0.166 


5 2- 

2 Z ' 


1+ 


78.60 


0.55(2) 


1.7(1) 


0.033 


- 1— 
2 1 ' 


2+ 


75.23 


0.43(2) 


1.8(1) 


0.033 


I 2 - 


2+ 


68.29 


0.30(2) 


1.8(1) 


0.299 


|2+, 


2- 


107.19 


0.95(2) 


1.5(1) 


0.132 


|3+, 


2- 


98.97 


0.89(2) 


1.5(1) 


0.009 


|2+, 


3- 


yo.oo 


U.oz^z ) 


1.6(1) 


u.uuy 


§3+, 


3- 


on r~\ 


U. / o[Z) 


1.6(1) 


U.loo 


- 3- 

2 J ' 


3+ 


112.38 


1.63(2) 


1.2(1) 


0.105 




4+ 


96.07 


1.56(2) 


1.2(1) 


0.136 


§4+, 


4- 


95.75 


3.22(4) 


0.36(7) 


0.086 


|5+, 


5- 


79.63 


3.45(4) 


0.23(7) 


0.105 


£5-, 


5+ 


52.81 


9.1(6) 


-2.6(3) 


0.072 




6+ 


36.85 


12.1(6) 


-4.1(3) 


0.085 


2 u ' 


6+ 


20.25 


-34.(2) 


19.(2) 


0.062 


13 7 
2 ' ' 


7+ 


36.06 


-18.(2) 


11.(2) 


0.071 


fr+, 


7- 


126.59 


-7.6(2) 


5.8(4) 


0.054 


¥ 8+, 


8- 


142.24 


-6.5(2) 


5.3(4) 


0.061 


¥»-, 


8+ 


268.76 


-4-7(1) 


4.4(3) 


0.047 


¥9-, 


9+ 


284.25 


-4-3(1) 


4.2(3) 


0.053 


f 9+, 


9- 


448.75 


-3.59(7) 


3.8(3) 


0.042 


f 10+, 


10- 


464.07 


-3.39(7) 


3.7(3) 


0.046 


¥ io- 


10+ 668.02 


-2.97(6) 


3.5(3) 


0.038 


2 ^ ' 


11+ 683.18 


-2.85(6) 


3.4(3) 


0.041 



There are three manifolds of levels, which belong to the 
vibronic states Il"x/2) n 3 / 2 , and S^, 2 (see Fig.jTJ). Accord- 
ing to (24 1 the strongest transitions take place between 



levels of the same manifold. The higher frequency transi- 
tions correspond to the change of the rotational quantum 
number J — » J + 1 (see Table IV I. Such transitions usu- 



ally have Q a » 
true for Z-C3H. 



0, Q A 



1 141 . We see that this is also 



For the II 1/2 and n 3 / 2 manifolds there is weak mono- 



A3/2 



tonic dependence of the sensitivities on J. This depen- 
dence is caused by the Coriolis interaction between these 
manifolds. For the upper part of the H3/2 spectrum we 
see some irregularities in sensitivities. They are caused 
by the resonant interactions with the nearby levels of the 
£1/2 manifold, where similar irregularities are observed 
for N > 22. All these irregularities are weak because 
interaction energy is much smaller than respective tran- 
sition frequencies. 

For the II states there are also lower frequency transi- 
tions between the levels of different parity with the same 



TABLE VI: Frequencies (MHz), Q 
(a.u.) for Tf-doubling transitions in 



factors and reduced MEs 
n 3 / 2 state. 



JF'p',Fp 


id 




Qa. 


<3m 


ii-Dir 


1 1-,1+ 


5 


61 


-2.63(8) 


3.2(2) 


1.493 


§2-1+ 


18 


50 


0.49(8) 


1.7(2) 


0.299 


§1-2+ 


-7 


30 


5.28(8) 


-0.6(2) 


0.299 


|2-,2+ 


5 


58 


-2.63(8) 


3.2(2) 


2.688 


|2+,2- 


22 


24 


-2.60(8) 


3.2(2) 


1.186 


§3+,2- 


31 


50 


-1.35(8) 


2.6(2) 


0.085 


|2+,3- 


12 


88 


-5.67(8) 


4.6(2) 


0.085 


§3+,3- 


22 


15 


-2.60(8) 


3.2(2) 


1.694 


f 3-3+ 


54 


92 


-2.57(8) 


3.2(2) 


0.943 


f 4-4+ 


54 


76 


-2.57(8) 


3.2(2) 


1.223 


9 4- 
2 + 


1 os 

lUo 


1 3 


^nfs i! 


0. 1 j 


1 o^n 


11 4- 

-2 — r 


1 

loO 


GO 

yy 


O A(\{R\ 
— Z.4D^o ) 


Q 1 (9\ 
0. 1 \Z) 


1 no7 
1 .uu / 


13 , 

T +~ 


9Q1 


71 


0/11 (Q\ 


^ 0^0^ 


OK/17 
U.o4 / 


15 , 

Tf — r 


/l 07 
4Z / 


Si7 


— Z.oO^o ) 


Q O^O^ 


O 707 
U. ( Z ( 


17 4- 

Tf +- 




Q 1 




n/o\ 
) 


O f^QQ 
U.OoO 


19 , 
T _ + 


70S 


ZO 




9 Ci(9\ 


tx p: 
U.OOo 


21 , 

T +~ 


1 HQ/1 
1U04 


ID 


O 01 {Q\ 

— z.zi ^y j 


9 Q(9\ 


O /1Q7 


23 , 

Tf r 


lOUO 


70 


O 1 7fQ >! 


9 $Z(9\ 


O AA« 


25 | 

Tf H 


1DUO 


07 

y 1 


O 1 ^(Q\ 

— z. lo^y ) 


9 ah ^ 


O /10^ 
U.4UO 


27 , 

T h 


iyoy 


Uo 


O 1 Q /^Q^ 


9 Q/1\ 

Z.oyl ) 


U.ODD 


29 , 

Tf +- 


oqoo 
zouu 


ID 


O 1 QfQ i! 

— z. lo^y ) 


9 7/11 


O QQ/1 
U.004 


31 , 

Tf r 


0fi8/1 
ZOo4 


y 1 




Z. 1 I 1J 


O ^Ofi 
U.OUO 


33 , 
Tf +~ 


OUoD 


ao 




Z. / yl) 


O 080 


35 

Tf — r 


o4yo 


R 1 




Z.O^l ) 


O Ofil 
U.ZD1 


37 , 
Tf + 


ooyo 


1 


^(~\ ^ 


z.ou^y ) 


O O/IO 
U.Z4Z 


39 1 

Tf — r 


4ZOO 


1 7 
1 I 


O Q/'1\ 

-z\y(lj 


Z.Oo^o ) 


O OO/I 
U.ZZ4 


- +- 
2 ~ 


4553 


04 


—3 5(1) 


2.42(5) 


0.208 


« , 
2 ~ 


4663 


43 


-4.6(2) 


2.2(1) 


0.192 


2 ~ 


4377 


16 


-7-5(2) 


1.4(3) 


0.174 


47 

2 T 


3097 


96 


-19.0(4) 


-2.3(9) 


0.149 


49 , 
2 "T" 


909 


06 


132.(2) 


53.(8) 


0.103 


61 , 
2 ~ 


19813 


69 


-3.11(5) 


-1.6(4) 


0.116 


53 , 
2 


16952 


67 


-1.31(2) 


0.0(4) 


0.136 


55 1 
2 "r 


16218 


56 


-0.61(2) 


0.8(4) 


0.138 



J (X-doublets). For diatomic radicals such transitions 
are known to be very sensitive to the variation of both 
constants [5HTU] . Electron spin gradually decouples from 
the molecular axis with growing rotational energy. As a 
result, the ^-doubling for low J values transforms to A- 
doubling for higher Js. In our case the electronic quan- 
tum number A is substituted by the vibronic quantum 
number K, otherwise the effects are rather similar (see 
Tables |v| and VI). Decoupling of the electron spin hap- 
pens around J = -y an d causes the anomaly in sensi- 
tivities for the II4/2 doublets around J = where the 
frequency drops below 50 MHz. For the Z-C3H molecule 
we can expect additional anomalies in sensitivities due to 
the proximity and strong interaction of II and S states 
One such anomaly is caused by the resonance be- 

Thc tran- 



tween the levels n 3 / 2 and E^y 2 with J 



49 
2 



sition frequency is higher here, about 1 GHz, but this is 
much smaller than for the neighboring rotational states. 

The hyperfine structure is much larger for the K- 
doublets of the H1/2 state. For this reason we do not 
neglect the hyperfine structure in Table [V] For high J 
values the transitions with AF ^ are strongly sup- 
pressed, so we list only transitions with AF = 0. In 



Table VI the hyperfine splitting is neglected for all but 
the first few transitions. For transitions with J > | the 
sensitivity coefficients for the hyperfine components of 
the transition are practically the same. 



TABLE VII: Frequencies (MHz), Q-factors and reduced MEs 



.u.) 


of 


some 


transitions IL./2 


Jp — > 


n 3 /2 J' P ■ 




./' 


P' 


./ 


P 




Q a 


Qm 


\\D\\ 2 


15 
2 




17 

2 


+ 


262072.96 


3.00(2^ 


-0.50(0) 


0.41 


17 

2 




17 

2 


+ 


456425.35 


1.70(2 


0.15(0) 


1.20 


1!) 

2 




17 

2 


+ 


674795.09 


1.14(1 


0.43(1) 


0.87 


15 
2 


+ 


17 

2 




261366.66 


3.02(2 


-0.51(0) 


0.40 


17 

2 


+ 


17 

2 




456743.26 


1.70(2 


0.15(0) 


1.20 


l!l 

2 


+ 


17 

2 




673718.38 


1.14(1 


0.43(0) 


0.84 


:il 

2 




33 
2 


+ 


185259.51 


3.45(2 


-0.71(1) 


2.31 


:',:! 
2 




33 

2 


+ 


558878.68 


1.14(1 


0.43(1) 


5.46 


35 
2 




:«:< 

2 


+ 


961254.71 


0.64(1 


0.68(1) 


3.30 


:il 
2 


+ 


:«:< 

2 




179719.83 


3.61(3' 


-0.81(1) 


2.24 


33 
2 


+ 


33 
2 




559110.84 


1.13(1 


0.43(1) 


5.49 


:',r, 
2 


+ 


33 
2 




954903.42 


0.66(i; 


0.67(0) 


3.20 


51 

2 




53 
2 


+ 


148111.71 


2.91(2 


-1.01(6) 


3.92 


53 
2 




53 
2 


+ 


729561.89 


0.67(1 


0.67(0) 


12.80 


55 
2 




53 
2 


+ 


1369461.89 


0.34(1 


0.82(1) 


6.63 


51 

2 


+ 


53 
2 




118586.34 


4.22(4 


-1.18(2) 


5.69 


.1:', 
2 


+ 


53 
2 




736802.86 


0.64(1 


0.63(1) 


10.93 


■IT) 
2 


+ 


53 

2 




1343531.64 


0.36(1 


0.81(1) 


7.01 



Because of the mixings (17) - (22 1 of the basic states 



there are also weaker transitions between n 3 / 2 , and 
SJ, 2 manifolds. Examples of such transitions are listed in 



Tables VII VIII Sensitivities of these transitions depend 
on the quantum numbers in a less regular manner, than 
sensitivities within each manifold. 

All transitions in Table |VII| have frequencies higher 
than 100 GHz. Because of that the sensitivity coefficients 
are not very high, but they are dispersed within intervals 
< Q a < 4 and -1 < < 1. Note that in order 
to study possible variations of fundamental constants we 
need to compare several transitions with different sensi- 
tivities. Thus, such a spread in sensitivities can be very 
useful HH. 



In Table VIII there are several low frequency transi- 
tions with very high sensitivities. Among them there are 
few with sufficiently high transition amplitudes. In par- 
ticular, there are three rather strong transitions at 27.6 
GHz, 25.1 GHz, and 10.5 GHz with sensitivities Q a from 
—5 to +19 and from —3 to +11. This is compara- 
ble to the sensitivities of the iVdoublet transitions from 
Tables |V| and |Vlj but for higher transition frequencies. 



TABLE VIII: Frequencies (MHz), Q-factors 
(a.u.) of some transitions H3/2 J p — > 
frequency means that final state lies lower. 



and reduced MEs 
NJ'p'. Negative 



N J' 


/ 

{] 


J 


P 






Qm 


\\D\r 


14 


2:1 
2 


4- 


31 
2 




— 159987 95 


2.00(2) 


1.90(7) 


0.10 


16 


31 

2 


4- 
— 


31 
2 




535601 40* 


—0.60(1) 


0.73(1) 


0.14 


16 


33 
2 


+ 


31 

2 





535512.07 


-0.60(1) 


0.73(1) 


0.11 


16 


33 
2 


+ 


33 
2 


- 


161892.90 


-1.96(2) 


0.12(7) 


0.21 


15 


:-»! 

2 


— 


53 
2 


+ 


-200166.07 


1.55(2) 


1.74(6) 


0.14 


17 


33 
2 


— 


33 
2 


+ 


540214.36* 


-0.58(1) 


0.73(2) 


0.18 


17 


35 
2 


— 


33 
2 


+ 


540229.47 


-0.57(1) 


0.73(2) 


0.15 


17 


35 
2 


— 


5.5 
2 


+ 


144436.89 


-2.13(2) 


-0.01(8) 


0.28 


16 


33 
2 


+ 


35 
2 


— 


-240483.13 


1.25(1) 


1.63(6) 


0.19 


18 


35 
2 


+ 


3o 
2 


— 


544660.22* 


-0.55(1) 


0.72(2) 


0.25 


18 


3 ( 
2 


+ 


3o 
2 


— 


544828.26 


-0.55(1) 


0.72(2) 


0.22 


23 


47 
2 




49 
2 


4 
— 


—514442 07 


0.17(0) 


1.20(2) 


7.13 


25 


49 
2 




■19 

2 


4 


579974 37 


—0.24(1) 


0.79(2) 


13.06 


25 


51 

2 




49 
2 


4 


569214.13* 


—0.15(0) 


0.90(1) 


6.18 


24 


47 
2 


4. 


49 
2 




18489 70 


— 13.9(2) 


-8 fl) 


0.15 


24 


49 
2 


4. 


■19 

2 




27624.21 


—5.28(7) 


—3.0(5) 


9.40 


26 


61 

2 


4. 


49 
2 




1162035.24 


—0.22(0) 


0.85(1) 


0.14 


25 


49 
2 




51 
2 


+ 


230.75 


-1099.(34) 


-742.(90) 


0.17 


25 


r>i 

2 




51 

2 


+ 


-10529.49 


18.8(2) 


11.(1) 


6.95 


27 


r,:>. 
2 




51 

2 


+ 


1188561.69 


-0.21(0) 


0.85(1) 


0.15 


24 


49 
2 


+ 


51 
2 




-571024.04* 


0.14(1) 


1.11(2) 


7.26 


26 


51 
2 


+ 


51 
2 




563386.99 


-0.34(1) 


0.79(3) 


10.93 


26 


5.", 
2 


+ 


51 

2 




556353.26 


-0.32(1) 


0.84(2) 


5.39 


26 


51 
2 


+ 


53 
2 




-18063.18 


13.6(2) 


11.(1) 


0.18 


28 


5:: 
2 


+ 


53 
2 




-25096.91 


9.2(1) 


6.7(8) 


3.01 


28 


55 
2 


+ 


53 
2 




1215046.58 


-0.20(0) 


0.86(2) 


0.17 



Transitions observed in Ref. |14] , 



Transitions observed in [14j are marked with asterisk in 
Tables |IV| and |VIII| All of them have frequencies above 
300 GHz and sensitivities, which are not very far from 
the typical rotational sensitivities: Q a ~ and Q M w 1. 
The maximal difference in sensitivities AQ a 0.7 and 
AQ^ « 0.4 corresponds to the transitions at 535.6 and 



571.0 GHz from Table |VIII| The only two transitions, 
which were detected at high redshifts in Ref. [30] are 
marked with the dagger in Table |IV| These transitions 
have much lower frequencies, but they are essentially ro- 
tational transitions with "normal" sensitivities. 

Let us discuss the accuracy of our calculations of the 
sensitivity coefficients Q a and Q^. As we mentioned 
above, we know the scalings of the parameters of the 
effective Hamiltonian with approximately 1% accuracy. 
So, we did several calculations of the sensitivity coeffi- 
cients. First, we changed all scalings by 1%. Second, 
we used scaling of the rotational constants from Eq. ( 29 ) 
keeping all other scalings unchanged. The uncertainty 
in the scalings of the smaller parameters of the effective 
Hamiltonian may be higher due to the nonadiabatic cor- 
rections. So we did two additional calculations with the 
scaling of the parameter /? changed by ±1/4, i.e. f3 ~ ^ 
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and P ~ fj, 3 / 2 . Finally, in order to check, how the fitting 
procedure may affect the results, we did several calcula- 
tions with slightly different sets of parameters. For ex- 
ample, we made a 13 parameter fit with three centrifugal 
corrections set to zero: Ju.d = Pd = Id = 0. In terms 
of the obtained frequencies, such fit is only three times 
less accurate than our final 16 parametric one. 

In Tables |IV| - |VIII| we give the average values of 
the Q factors for all calculations, described above. 
The errors, given in the brackets, correspond to the 
maximum deviations from these average values for 
individual calculations. In most cases these errors are 
smaller than, or of the order of 10%, even for the large 
sensitivities. This accuracy is sufficient for the analysis 
of the experimental and observational data. 



sionless ratio £ — k' /k of the force constants in the two 
perpendicular planes, which include molecular axis. Pa- 
rameter E does not depend on the fundamental constants 
and vibrational intervals scale in the same way as for har- 
monic oscillator, i.e. E v ~ /i 1 / 2 . However, the Renner- 
Teller interaction modifies vibrational spectrum and can 
lead to the close lying vibrational states. Such states 
then strongly interact with each other due to the Cori- 
olis interaction. As a result, the molecules with strong 
Renner- Teller interaction can have low frequency mixed 
ro-vibronic transitions with strongly enhanced sensitivity 
coefficients to the variation of a and /i. For the Z-C3H 
molecule we found several types of transitions with sen- 
sitivity coefficients varying in a wide range. This opens 
new possibilities to study variation of fundamental con- 
stants in astrophysics. 



VI. CONCLUSION 

We have studied sensitivity coefficients to the varia- 
tion of the fundamental constants a and fx for the mi- 
crowave and submillimeter spectra of the linear poly- 
atomic molecule with strong Renner- Teller interaction. 
As an example we have chosen Z-C3H molecule, which is 
often observed in the interstellar molecular clouds and 
recently has been detected at the redshift z — 0.89. 

The Renner- Teller interaction depends on the dimen- 
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